Linear Responses in Time-dependent Hartree-Fock-Bogoliubov Method with Gogny 
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A numerical method to integrate the time-dependent Hartree-Fock Bogoliubov (TDHFB) equa- 
tions with Gogny interaction is proposed. The feasibility of the TDHFB code is illustrated by the 
conservation of the energy, particle numbers, and center-of-mass in the small amplitude vibrations 
of ^°0. The TDHFB code is applied to the isoscalar quadrupole and/or isovector dipole vibrations 
in the linear (small amplitude) region in oxygen isotopes 18'20,22,24q^ titanium isotopes **4,50,52,54rpj^ 
neon isotope Ne, and magnesium isotopes ' Mg. The isoscalar quadrupole and isovector dipole 
strength functions are calculated from the expectation values of the isoscalar quadrupole and isovec- 
tor dipole moments. 

PACS numbers: 21.60.-n 



I. INTRODUCTION 

In nuclear physics, the mean-field theory has played 
a central role in understanding the nuclear ground state 
properties as well as the dynamical nuclear excitations 
not only with small amplitude but larger one. Amomg 
the microscopic methods in the nuclear mean-field theory, 
the time-dependent Hartree-Fock (TDHF) is character- 
ized by the full self-consistency between the mean-field 
(collective) motion and the internal nucleonic motion ill] . 

The first practical numerical calculation of the TDHF 
was carried out more than thirty years ago [2]. Since 
that time, the TDHF has been applied to the heavy- 
ion reaction processes such as fusion, deep-inelastic colli- 
sion, fission and so on with the purpose of understanding 
the microscopic mechanism of the large amplitude nu- 
clear collective motions in the collision processes [^-Q- 
In the TDHF calculations, the three-dimensional (3D) 
Cartesian grid has been used, combined mainly with the 
Skyrme effective interactions. Recent TDHF calculations 
of the reactions of the heavier nuclei are carried out by 
taking advantage of the rapid increase of the computer 
power [i4i|. 

The TDHF has been widely used also in the small am- 
plitude domain of the collective motions such as giant 
resonances (GR) . It is well known that the TDHF is con- 
nected with the random phase approximation (RPA) in 
the small amplitude limit, which is a useful tool to study 
the collective excitations of the nucleus. Combined with 
the linear response theory, the excitation energies of the 
collective motions are obtained by calculating the small 
amplitude responses of the nucleus to the impulse type 
perturbation as an initial condition of the TDHF equa- 
tions Q- A new method to calculate the RPA energies 
and am plit udes was proposed and applied to several nu- 
clei [ISlnl. 

The quasiparticle RPA (QRPA) is an extention of 
the RPA to include the pairing correlation effects. In 
this decade, with the purpose of studying the excitation 



mechanism of the unstable neutron rich nuclei, the QRPA 
calculations have been carried out based upon the various 
formulations ^^^. 

Just as in the case of the RPA as a small amplitude 
limit of the TDHF, the QRPA is the small amplitude 
limit of the time-dependent Hartree-Fock Bogoliubov 
(TDHFB). The TDHFB is an extension of the TDHF, 
in which the pairing correlation in the nucleus is treated 
selfconsistently with the particle-hole (ph) correlations 
which are taken into account in the TDHF. 

A possibility of the practical numerical calculations 
of the TDHFB with an effective nuclear force (Gogny 
force) was presented in Ref. [28||. The type of full TD- 
HFB calculation was carried out in the study of the pair- 
ing vibrations with Skyrme interaction [29|. Recently, 
a framework of the TDHFB calculation in the canonical 
form with the Skyrme interaction (Cb-TDHFB) was pro- 
posed and applied to the linear responses of both light 
and heavy nuclei [S^- The cpu time is substantially re- 
duced even in the calculations of the heavier nuclei by 
using the Cb-TDHFB, which is a useful method in study- 
ing the systematic property of the nulei in the wide range 
of the nuclear chart. 

In the HFB calculations with the Skyrme interactions 
(Skyrme-HFB), the Skyrme interaction is used for the ph 
channel, while another type of interaction is introduced 
for the particle-particle (pp) channel. Since the zero- 
range interaction is assumed in the Skyrme HFB, it is 
necessary to set the appropriate cut-off energy and choose 
the optimum parameter set in the pairing part [Slj . 

In contrast with the Skyrme HFB, in the HFB cal- 
culations with the Gogny interaction (Gogny- HFB), the 
ph channel and the pp channel are treated on an equal 
footing. In the Gogny-(TD)HFB, a practical cut-off of 
the energy range of the physical space is introduced nat- 
urally from the terms with the Gaussian functions of the 
Gogny interaction. Therefore, the Gogny interaction is 
suitable for the formulation of the self-consistent TDHFB 
framework combined with a practical numerical method 
of integrating the TDHFB equations. We put emphasis 



on this advantageous point of the Gogny interaction and 
developed a Gogny-TDHFB code. 

In this article, we would like to report the first results 
of the Gogny-TDHFB calculations, illustrating the fea- 
sibility of the formulation and numerical codes. In sect. 
2, the basic equations of the TDHFB are simply given. 
A method of solving the TDHFB equations is described. 
In sect. 3, taking oxygen isotope ^'^O as an example, 
the feasibility of the numerical code of integrating the 
TDHFB equations is illustrated. The conservations of 
total energy, nucleon numbers, and center-of-mass are 
discussed. In sect. 4, the Gogny-TDHFB is applied to 
the linear (small- amplitude) vibrations of spherical nuclei 

18,20,22,240^ 26Ng^ ^ud 44,50,52,54rpi ^g ^^^ ^^ deformed 

nuclei ^"'''^'^Mg as typical examples. The strength func- 
tions of the isoscalar quadrupole and isovector dipole ex- 
citations are calculated from the linear responses. Sec- 
tion 5 is for summary and concluding remarks. 



II. TDHFB EQUATION 

The nuclear Hamiltonian is assumed to take the form, 
H = 22 TapC'lCp + T 2_/ '^aMsCl^Cf^CsC^, (1) 



Q/3 



a/97(5 



where Tajj is kinetic energy matrix element and Vajj-yS is 
antisymmetrized two-body matrix element. The operator 
C^(Cq) is a nucleon creation (annihilation) operator of 
a state labelled with a. 

The TDHFB equation for the generalized density ma- 
trix TZ is written in the form [l|, 



ihn= [H,n], 

with the HFB Hamiltonian T-L, 



n = 



h A 



-A* -h* ' ' 
and the generalized density matrix 7^, 



-JZ 



-K* 1- p* 



(2) 



(3) 



(4) 



where p and k are normal density matrix and pairing 
tensor, respectively. 



Pc.0 = {V*V^) ^„^ = (V^*[/^) 



af3 ■ 



(5) 



The matrices U and V are the Bogoliubov transformation 
matrices. 



a 



(6) 
(7) 



from the particle operators C^ and Cq into the quasi- 
particles /3jJ and Pk- The notation U^ {V'^) stands for the 
transposed matrix of U (V). The mean field Hamiltonian 
h and the pairing mean field A are introduced through 
the relations. 



h 



afi 



T, 



ap 



Q/3, 



(8) 

^ap = ^^^a^pSPS'yi A^^ = - 2^ Va/J-y^K^a, (9) 
7(5 

with the normal density p and the pairing tensor k. 

Instead of the equation ^ of the generalized den- 
sity matrix 7^, we write the equation of motion in the 
form [H,!!! 



,n^ ( uit) \_^( u{t) \ 
'^dt \ v{t) )-^\ V{t) ) ' 



(10) 



for the matrices U and V in the Bogoliubov transforma- 
tion ® and 0. 

Noting that the form of the equation ^TO\i is analo- 
gous to the TDHF equation [33| , we follow the numerical 
method of integrating the TDHF equation and get the 
formal solution of the TDHFB equation (ITI 



(«+i) 



exp 



_,^-^(n+l/2) 

h 



(") 



(11) 



where 



(") 



is a matrix composed of matrices U and 



y at a discretized time t„ {n = 0, 1, 2, • • • ). The quan- 
tity 7^("+i/2) is an adequate TDHFB Hamiltonian ma- 
trix at every time step from t„ to i„+i = tn + At. The 
Hamiltonian matrix 7^("+i/2) js determined so that the 
expectation value E of the Hamiltonian ([1]) with respect 
to the time-dependent HFB state |<I>), 



($|i/|$) 

a/3 



1 

2^^ 



r 



a0^al3, 



fl2) 



is kept constant in the course of the numerical integra- 
tions in ([TT|). 

Since in the TDHFB calculation the energy conserva- 
tion is one of the most important conditions to be ful- 
filled, let us see how the energy in (fT2|) is conserved with 
respect to a small variation in the matrices U and V. 

Let us assume the time step At is small enough to 
approximate the right hand side of Eq. (fTTT) within the 
first order of At in the power series expansion of the 
exponential function. Putting a parameter A = ■^, the 
variations of the matrices U and V at a time f into the 
new matrices U' and V' are written within the first order 
of the parameter A, 



U' 
V 



-iX 



h 

-A* 



A 



U -i\{hU + AV) 
V + iX{h*V + A*U) 



(13) 



Here, for the ease of the discussion, let us assume also 
that the time increment Ai is small enough so that we can 
identify the mean field Hamiltonian /i("+i/2) and mean 
pairing potential A("+^/^) in the TDHFB Hamiltonian 
^(n+i/2) ^[fj^ ]j g^jj^fj /^ g^^ ^j^g time t, respectively. 

Using the relations in Eq. (|13p . the variations in the 
density matrix p and pairing tensor k in Eq. ([S]) are 
expressed, respectively, as 

p' = V'*V'^ 
= {V* ~ iX (AC/* + hV*)) {V^ + iX (t/^At + V^h^)) 
= p-iX[h,p]~iX{-AK* + K.A*), (14) 

k' = V'*U'^ 
= (V* - tX (At/* + hV*)) {U^ - ^A {U^h'^ + F^A^)) 
= K-iX{AU*U^ + hn + Kh* - pA) . (15) 

We see that the conservation of the particle number in 
every time step holds from Eq. (|14p : Taking the trace of 
both sides in Eq. ([T4)) . we have 



Tip' — Tip 

-iX (Tr {[h, p]) + Tr {-Ak* + kA*)) 



i-Z 2_^ Z_^ [yal3fj.u'ifii''^al3 I^a0i^ij.va0l^^iv) 



afS fJ.1^ 



= 0, 



(16) 



where the definition of A in Eq. ^ and the anti- 
symmetric property of the matrix elements Vq.^^^ and 
pairing tensor Kap are used. The notation Tr stands for 
taking the trace of the matrices. Then, the particle num- 
ber Nr = Tip is conserved from step to step. Here, Nr is 
the proton (neutron) number Z (N) for proton (neutron) 
density matrix p — Pp {pn), respectively. 

Putting the expressions of the density and pairing ten- 
sor up to the first order in the parameter A in Eqs. (|14p 
and (fT5)) into the expression of the energy in Eq. ([T^. 
we have the variation of the energy as follows: 



SE = TT{~iXh[h,p] ~iXh{-AK* + kA*)} 

~Y 

_iX^ 

~ 2 

= 0, (17) 



-Tr{-A*A - A*hK - A*Kh* + A*Ap* + A*pA} 
-Tr{ A*A + h*K*A + K*hA~- A*pA- p*A*A} 



where the relations h* — h^ , k^ ~ — k, and A-^ = —A 
are used. From Eq. dTT]) . we see that we can inte- 
grate the TDHFB equation ()10p . conserving the energy 
expressed as in Eq. ()12p . by appropriately setting the 
time increment At and intermediate TDHFB Hamilo- 
nian H^""*"^/^-*. In the case of the Gogny interaction 
[SJ, |35| with the parameter set DIS, the contribution of 
the density-dependent term is included only in the mean- 
field Hamiltonian h in the HFB Hamiltonian ([3]) , through 



the variation of the nucleon density p(x) [3^. Then, the 
energy conservation relation p7)) holds when the Gogny 
interaction with the parameter set DIS is adopted in the 
Hamiltonian ([T|). 



III. NUMERICAL SOLUTION AND 
CONSERVED QUANTITIES 

A. Parameters 

With the purpose of demonstrating the feasibility of 
our TDHFB code, we apply them to the small amplitude 
(linear) responses of the oxygen isotope ^°0. We adopt 
the Gogny interaction [3J, [331 ^O"^ the two-particle inter- 
action part in the Hamiltonian ([T]). The parameter set 
DIS is used. The Coulomb part is omitted in all of the 
calculations in this article. 

The initial condition in the present calculation is of the 
impulse type: The static HFB solutions Uq and Vb are 
changed into the initial matrices C/^^-' and y(°) by the 
relations 



N„ 



l/(°) = exp(z£Q)Vb= Y^ 
C/(o) = cxp(-^eQ*)f/o= ^ 



i'^e'^Q 



■Vo, 



(18) 



i^i-eyq*" 



-Uo, 



(19) 



where Q stands for the matrix representation of a multi- 
pole operator with respect to the numerical basis states. 
Nmax is the maximum number up to which the exponen- 
tial function is expanded into power series. The param- 
eter e is chosen so that the linearity is held in the course 
of the time integration of the equations (TTOl) . 
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FIG. 1: Time dependence of expectation value of isoscalar 
quadrupole moment of O. 

In making the TDHFB hamiltonian H, we used the 
same form of energy functional as that used in the calcu- 
lation of the HFB ground state, i.e., the center-of-mass 
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FIG. 2: Deviation of the excitation energy from the initial (t 
= 0) value of the isoscalar quadrupole vibration in fig. [l]with 
respect to time. 
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FIG. 3: Deviation of neutron (proton) number expecta- 
tion value from the accurate number 12 (8), respectively, in 
isoscalar quadrupole vibration in fig. [T] 



correction part Ecm ~ —P'^/2Am witli tlie nucleon mass 
m and tlie total momentum operator P = X^i Pi oi A 
nucleons is included. 

As the numerical basis, the three-dimensional har- 
monic oscillator eigenstates were used. The two-particle 
matrix elements Va^-yS in the Hamiltonian ^ were cal- 
culated by following the method described by Girod and 
Grammaticos j36i] . 

We used the space of basis states with the relation 
< Tix + Uy + Uz < 4, where rix {uy, Uz) is the number 
of quanta of the harmonic oscillator basis states in the 
X (y, z) direction, respectively. The harmonic oscillator 
parameters were set as follows: fiiiJx = huoy = fkOz = 13.7 
MeV. 

In the initial conditions ([18]) and ([T9| , we took the mul- 
tipole operator Q to be an isoscalar quadrupole operator 
Qap — (2z^ — x"^ — y^)^g- The parameter e was put to 
be 1.0 X 10"^, being small enough so that the linearity of 
the oscillation with respect to e is satisfied. The power 
series expansions of the exponential functions were taken 



up to the Nmax-th. order with N^ax put to be ten. 

In the course of the time integrations in pTj). a 
predictor-corrector method was made use of in making 
the Hamiltonian matrix 7{("+i/2) at every time step [331. 



B. Conservation of energy, particle number, and 
center-of-mass 



In fig. [TJ we display the time variation of the isoscalar 
quadrupole moment Q2o{t) = {J2i=i'^iY2o{fi)) of ^"O. 
The time increment cAt is 0.2 fm. After the initial im- 
pulse, regular small-amplitude oscillations take place. 

In fig. [21 the deviation of the excitation energy from 
the initial (t = 0) value in the course of the vibration 
is shown. The excitation energy varies from the initial 
value 0.02554 MeV to 0.02561 MeV after integration time 
4500 fm (22500 steps). The deviation of the excitation 
energy from the initial value in the course of the time 
integration is reduced when we carry out more numbers 
of predictor-corrector loops in each time step. 

In fig. [21 the deviation of the expectation values of the 
nucleon number Nr from the accurate values (8 protons 
and 12 neutrons) is displayed with respect to time. In 
^°0, the protons are in the normal state, whereas the 
pairing correlation is active on the neutron side. There- 
fore, the integration of the equations of motion of the 
proton orbitals is equivalent to the TDHF case. Then, 
the total number of protons are conserved within 10"^^. 
The neutron number in fig. [31 is kept within around 10~^ 
in the present integration process. These results illustrate 
that the unitarity of the time-displacement operator in 
Eq. (ITTI) holds within this accuracy. 




Pi 



s o /|^^;Vl^fVvVW^ 



Time (fm) 



FIG. 4: Time dependence of expectation values of isovector 
dipole moment D(t) (upper panel) and center-of-mass R(t) 



(lower panel). 



Taking the initial impulse operator Qap to be the 
isovector dipole operator Dap — ^ {^)a0 ^'-"' P^'otons 
and —4 (•z)q,« for neutrons, we have the oscillation pat- 
tern in fig. [31 where the expectation value D{t) = 



■E- 



■ J2i=i ^i) is displayed together with the experimental results. 



center-of-mass expectation value in z direction R{t) 

^{J2i=i^i + J2i=i^i) i'^ the first 12000 steps of inte- 
gration. Here, N{Z) stands for the neutron (proton) 
number, respectively. We see that the fluctuation of the 
center-of-mass motion R{t) is kept within 2 x 10^^ fm. 
The effect of the center-of-mass motion on the isovector 
dipole moment D{t) could be neglected. 

When the quadrupole moment operator (2z^ — x^ — ij^ 
or x^ — y^) is used as the impulse operator, the expecta- 
tion values of the center-of-mass and off-diagonal terms 
xy, yz, and zx are kept to be zero, since the symmetry 
of the system at the initial time i = is kept throughout 
the integration time. 



IV. APPLICATIONS 



A. Spherical nuclei 



18,20,22,24 



o 




10 20 30 10 20 30 40 



Energy (MeV) 



Using the impulse type initial conditions (Eqs. ([18]) 
and (J19p ) with an isoscalar quadrupole operator 
Qap — (2z^ — x"^ ~ y'^) aB ^^'^ ^^^ space of basis func- 
tions Q < rix -\- n^ + TLz < 4, we applied the method of 



"y 



integrating the TDHFB equations discussed in the pre- 
vious section to the oxygen isotopes i8.20,22,24q_ ^q^- 
ing the Fourier transformation of the quadrupole mo- 
ment Q2o{t), we calculated the strength functions of the 
isoscalar quadrupole mode of the oxygen isotopes in fig. 
[5]@. 

In fig. [SI we see that the low-energy peaks are located 
below 5 MeV. The energies of the peaks below 5 MeV are 
shown in fig. [Hi The neutron pairing energies of the oxy- 
gen isotopes are shown in table 1. In the four oxygen iso- 
topes 18.20,22,240^ |-]-^g magnitude of the neutron pairing 
energy is the largest in ^"O and zero in ^''O. The neutron 
number 16 in ^''O corresponds to the subshell closure of 
2s| . The tendency of the energies of the peaks in the TD- 
HFB calculation that the lowest (highest) is of '^°0_^'^0) 
is in accordance with the experimental data [37(-[39[. In 
the present TDHFB calculation, the energy of the lowest 
peak of ^*0 is 4.1 MeV, which is similar to the QRPA 
predictions around 4 MeV [l^ Sol- The experimental 
results of the energies of the lowest 2+ state in ^^O are 
around 4.7 MeV [sl,!!^, which is higher than the QRPA 
and present TDHFB results by around 0.7 MeV. 

Because the Coulomb force is omitted in the present 
calculations, we compare our result of the first 2+ energy 
of the oxygen ^®0 with the experimental value of neon 
^^Ne. The experimental values of the first 2+ energies of 
the oxygen ^®0 and neon ^^Ne are 1.98 MeV and 1.89 
MeV 4l|, respectively. These two nuclei are mirror nu- 
clei to each other and their first 2+ energies are located 
within the width of 0.1 MeV. The present TDHFB resuh 
of the oxygen ^^O is around 0.5 MeV higher than these 



FIG. 5: Strength functions of isoscalar quadrupole mode of 

1 A 0C\ 00 A. 

oxygen isotopes ' ' ' O. Artificial width of 1 MeV is used. 

Through the Fourier transformation of the isovector 
dipole moment D{t) of ^°0 in fig. [4| we get the strength 
function of the isovector dipole mode in fig. [7| The peak 
is located at 25 MeV, which is larger than the observed 
value by a few MeV [42] . Since the absolute values of the 
center-of-mass is kept practically small in the time inte- 
gration, the effect of the spurious motion in the strength 
function might be neglected. 
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FIG. 6: Energies of the first peaks of the isoscalar quadrupole 
strength functions of oxygen isotopes i8,20,22,24q -^^ ^^ |gj 
Closed (open) squares stand for the TDHFB (experimental) 
results, respectively. Plus (-(-) is the first 2^ energy of neon 



TABLE I: Neutron pairing energy En in oxygen isotopes 

16,18,20,22,24q -^ Gogny_HFB. 



En (MeV) 



^O i«o 



^O 



2q 24q 



- 4.56 - 5.33 - 2.39 




20 30 

Energy (MeV) 



FIG. 7: Strength function of isovector dipole mode of ' 
Artificial width is 0.5MeV. 



^O. 



harmonic oscillator eigenfunctions. We need some theo- 
retical tools to separate out the effects of the pairing cor- 
relation. We expect that the idea of the Cb-TDHFB will 
be a good starting point to improve the present formula- 
tion and attack the problem of making clear the pairing 
effects in the dynamical processes of nuclear collective vi- 
brations not only of the quadrupole type but of the other 
multipolarity types such as monopole and octupole. 
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2. Pairing energy in O 



In fig. m we show the time dependence of the pair- 



ing energy Ep 



(Gauss) 



4 A^a 



». , .(Gauss) 



K^y 111 of the 



Gaussian part V)^J™^ in the Gogny interaction in the 
quadrupole oscillations in fig. [T]of ^°0. Here, k.^„ is the 
pairing tensor of the neutrons. The magnitude of the 
pairing energy which comes from the other parts in the 
Gogny interaction, i.e., spin-orbit, center-of-mass, and 
Coulomb (if included in the calculations) terms, is small 
compared with that from the Gaussian part. Since the 
initial impulse is taken to be small enough to ensure the 
linearity of the quadrupole oscillation, the fluctuation of 
the pairing energy in flg. |8] is within 0.05 MeV in the 
integration time up to 4500 fm. 

We take the Fourier transformation of the fluctuation 
of the pairing energy and show the Fourier (sine) ampli- 
tudes in the lower panel in fig. [5] Corresponding to the 
giant resonance peak in the strength function in the en- 
ergy range from 20 MeV to 25 MeV, we see remarkable 
fluctuations in the Fourier amplitudes. In the low energy 
region under 7 MeV, we flnd fluctuations in the Fourier 
amplitudes corrresponding to the low-energy peak at 2.3 
MeV and a shoulder around 4.4 MeV in the strength 
function. On the other hand, there are fluctuations in the 
Fourier amplitudes to which there are no corresponding 
peaks in the strength function curve. 

In the present formulation of the TDHFB, we can not 
separate out the effects of the pairing correlation, since 
they are embedded in the dynamical changes of the ma- 
trices U and V which are deflned on the fixed basis of 3D 



FIG. 8: (Upper panel) Time dependence of the pairing en- 
ergy Ep ^^ of the Gaussian part in the Gogny interaction 
in the quadrupole oscillation in fig. [T]of ^"O. (Lower panel) 
Fourier (sine) amplitude of the fluctuations of the pairing en- 



ergy E), 



(Ga 



s) 



in the upper panel (solid curve, left axis). Ar- 



tificial width of 0.25 MeV is used. The isoscalar quadrupole 
strength function is plotted in dashed curve (right axis). 



q 44,50,52,54 jij 

The next example of the TDHFB calculations in the 
spherical nuclei is the application of the method to the 
isoscalar quadrupole vibrations in the titanium isotopes 
44,50,52,54rpj^ Here and hereafter, the space of the basis 
functions < n^. + rij, + n^ < 5 was used in the calcula- 
tions. 

The strength functions of the titanium isotopes are 
shown in fig. 1^1 From the strength functions, we got the 
energies of the first low-energy peaks, which are shown 
in fig. [TUl In fig. [ini we see that the TDHFB results of 
the lowest-energy 2^ peaks are along the experimental 
values [il,!!!]: The lowest 2+ energies in 5°Ti and ^-^Ti 
are higher than the other two nuclei '*'*Ti and ^^Ti. 

This tendency of the lowest 2+ energies of the four nu- 
clei is in accordance with the tendency of the neutron 
pairing energies and sushell closures. In table 2, we show 
the pairing energies of protons and neutrons in the four 
nuclei 44,50,52,54rpj jj^ ^j^gjj. jjpg gj.oyj^(i states. The pro- 
ton pairing energies in 44,50,52,54rpj ^-^^ almost constant 
around - 4.8 MeV. In ^°Ti, the neutron number is 28 and 



lf| orbitals are filled. The pairing energy in ^'^Ti is zero. 
Likely, in ^''Ti, 2p| orbitals are filled and the pairing 
energy is comparatively small. 




10 20 30 10 20 30 40 



Energy (MeV) 



TABLE IL Proton (neutron) pairing energy Ep (En) in tita- 
nium isotopes ■44.50,52,54rp. j^ Qogny-HEB. 



Ep (MeV) 
En (MeV) 



'Ti ^''Ti 



■"Ti 



*Ti 



4.92 - 4.72 - 4.79 - 4.80 
4.92 - 2.75 - L17 



^^Ne is one of the neutron-rich nuclei, and a pygmy state 
is known as the low-energy components in the isovector 
dipole resonance |45| . 

Taking the isovector dipole operator as the impulse op- 
erator, we calculated the expectation values of the isovec- 
tor dipole moment in ^^Ne. We show the strength func- 
tion in fig. [TTJ 

The energy of the small 1^ peak near 10 MeV is 11.9 
MeV in the present TDHFB calculation. In the QRPA 
calculation by Peru at al. [22|, they got the value 10.64 
MeV for the 1~ peak of the pygmy state. The experi- 
mental value of the pygmy state is around 9 MeV [45| . 
The corresponding energy by the present TDHFB calcu- 
lation is larger than the QRPA (experimental) result by 
1.3 (3) MeV, respectively. 



FIG. 9: Strength functions of isoscalar quadrupole mode of 
titanium isotopes 44,50,52,54r[.j Artificial width of 1 MeV is 
used. 
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FIG. 11: Strength function of isovector dipole mode of ^^Ne. 



B. Deformed nuclei 



FIG. 10: Energies of the first peaks of the strength functions 
of the isoscalar quadrupole mode of the titanium isotopes 
44,50, 52,54rpj j^ g^^ [g] Qj^gg^j (open) squares stand for the 

TDHFB (experimental) results, respectively [il. li^. 



4- 



'Ne 



The mean-field of the ground state of ^^Ne is spherical 
in the Gogny-HFB calclation. The protons are in the 
paired state whereas the neutron pairing energy vanishes. 



As the typical examples of the deformed nuclei, we 
took magnesium isotopes ^'^Mg and '^'*Mg. In ^'^Mg, the 
pairing energies of both protons and neutrons are zero 
in the HFB ground state, in which the mean-field is in 
a prolate form. ^''Mg is also prolately deformed in the 
HFB ground state, in which the neutron pairing energy 
is finite. The deformation parameter /3 of ^'^Mg ('^'*Mg) 
is /3 = 0.51 (0.43), respectively 0]. 

Calculating the time dependence of the expectation 
values of the K = and 2 components of the isoscalar 
quadrupole moment, we get the strength functions in fig. 
[T^ Here we used the isoscalar quadrupole operator x'^ — 



y^ as the initial impulse operator to excite the K ~ 2 though the former tend to be smaller than the latter. 



component of the isoscalar quadrupole vibrations. 

The giant resonance peaks with K = ^ and 2 are seen 
in the energy region from 15 MeV to 30 MeV in both of 
the nuclei ^^Mg and ^-^Mg. The peaks with X = 2 are 
located in the higher energy region than those with K = 
0, reflecting the prolate deformation of the mean- fields in 
the ground states of the two nuclei. 

The energy of the lowest K = (2) peak in ^^Mg is 
8.06 (3.3) MeV, which is comparable with an experimen- 
tal value 6.43 (4.23) MeV, respectively [13 . The charac- 
teristic point that the K — 2 peak comes lower in energy 
than that with i^ = is closely related with the soft- 
ness of ^''^Mg with respect to the direction of the gamma 
deformation. 

In ■^'*Mg, the energy of the lowest K = (2) peak 
is 2.07 (3.0) MeV. In the QRPA calculations with the 
Skyrme interaction, the energy of the corresponding peak 
is 2.65 (around 3.0) MeV 24, 27]. The energies of the 
two peaks with K — and 2 by the TDHFB calculations 
are comparable with those by the QRPA calculations. 



The strength functions of the isovector dipole modes 
of ^^Mg are shown in the left panel in fig. [T31 The en- 
ergy of the a: = (1) peak in fig. [His 20.2 (28.9) MeV, 
respectively. The corresponding value obtained by using 
the QRPA by Peru et al. is 19.35 (27.44) MeV [25] . The 
present TDHFB result of the K = {1) peak energy is 
higher by 0.9 (1.5) MeV than the QRPA result, respec- 
tively. 

The strength functions of the isovector dipole modes 
of ^^Mg are shown in the left panel in fig. [T31 The en- 
ergy of the Js: = (1) peak in fig. [13] is 20.2 (28.9) MeV, 
respectively. The corresponding value obtained by using 
the QRPA by Peru et al. is 19.35 (27.44) MeV [H]. The 
present TDHFB result of the K ^ (1) peak energy is 
higher by 0.9 (1.5) MeV than the QRPA resuh, respec- 
tively. 

From the isovector dipole strength functions, we cal- 
culated the photoabsorption cross section a{E) based on 
the strong coupling scheme [3, ]^ 



J 



aiE) 
dE 



47rV^ 

ch 
1 

irhe 



K=-l 



dBjEJiK) 
dE 



-Im 



dt{^{t)\hK\^it))e 



iEt/h--it 



(20) 



(21) 



where 



dBjEJiK) 
dE 



is a transition strength with isovector V. SUMMARY AND CONCLUDING REMARKS 



T J2i=i ii^YiK)i for protons and 



for neutrons with K = ±1,0. 



dipole operators fix 

The nuclear state |*(t)) is the TDHFB state. 7 is an 
artificial width used in the strength function and e is the 
small parameter used in the initial conditions (|18|) and 
(dSl) of the TDHFB equations. 



In the right panel in fig. [T3]we compare the photoab- 
sorption cross section with the experimental results. The 
K = peak comes near the one of the experimental cross 
section, whereas the jiiT) = 1 peak is located in the higher 
energy region than the distribution of the peaks in the 
experimental data [49|]. 

The difference of the TDHFB results presented above 
in this section from the RPA results and experimental 
data might come partly from the restricted space of the 
basis functions {rix + Uy + Uz < 5) and the omission of 
the Coulomb force in the present TDHFB calculations. 
It is necessary to study the effects of the size of the space 
of basis functions and Coulomb force on the TDHFB re- 
sults. 



In this article, we presented a practical numerical 
method of integrating the TDHFB equations with Gogny 
interaction. The three-dimensional harmonic oscillator 
eigen functions were used as the numerical basis func- 
tions. The conservation property of the energy, parti- 
cle numbers, and center-of-mass of the method was il- 
lustrated. The method was applied to several typical 
spherical nuclei as well as axially deformed ones and the 
strength functions were obtained in the framework of the 
linear response theory. It was shown that the tendency 
of the energies of the low-energy peaks in the present 
TDHFB calculations was similar to the one of the ex- 
perimental results, though the energies obtained in the 
TDHFB calculations tend to be lower than the QRPA 
and experimental results. On the other hand, the en- 
ergies of the typical peaks in the energy region of the 
giant resonances tend to be larger than the QRPA and 
experimental results. We will study carefully how the 
numerical results depend on the space of basis functions 
and Coulomb force. 

The present calculations were carried out under the 
conditions: i) The maximum major shell quantum num- 
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FIG. 12: Strength functions of isoscalar quadrupole mode of 
magnesium isotopes ■^*'^''Mg. Solid (dashed) curve stands for 
K = {2) component, respectively. Artificial width of 1 MeV 
is used. 
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FIG. 13: Strength functions of isovector dipole mode (a)) and 
photoabsorption cross sections (b)) of ^^Mg. Solid (dashed) 
curve stands for \K\ — (1) component, respectively. Artifi- 
cial width of 1 MeV is used in the strength functions. Exper- 
imental data of cross sections are shown with plus (-I-) [49|]. 



ber of the space of the harmonic oscillator basis func- 
tions was restricted to be four or five, ii) The Coulomb 
force was omitted in both HFB and TDHFB calculations. 
These restrictions were introduced to save the cpu time 
and check the feasibility of the present numerical codes 
of TDHFB. 



It is necessary to use larger space of basis functions 
and to include the Coulomb force so that we can deal 
with the problems of heavier, unstable nuclei in terms of 
the present TDHFB method. When the larger space of 
the basis functions is taken, we need efficient numerical 
methods to save the cpu time of the iterations in the 
full Gogny-TDHFB calculations with the Coulomb force. 
The acceleration of the calculations of the Coulomb part 
might be realized by introducing a set of small number of 
optimal integration points and weights in the calculations 
of the Coulomb mean-field matrix elements. The new 
numerical codes by which we could save the cpu time are 
now being developed. 
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